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We report observation of crystallization of the glass- forming binary Lennard-Jones liquid first used 
by Wahnstrom [G. Wahnstrom, Phys. Rev. A 44, 3752 (1991)]. Molecular dynamics simulations 
of the met ast able liquid on a timescale of microseconds were performed. The liquid crystallized 
spontaneously. The crystal structure was identified as MgZn2. Formation of transient crystallites is 
observed in the liquid. The crystallization is investigate at different temperatures and compositions. 
At high temperature the rate of crystallite formation is the limiting factor, while at low temperature 
the limiting factor is growth rate. The melting temperature of the crystal is estimated to be 
Tm = 0.93 at p = 0.82. The maximum crystallization rate of the A2B composition is T = 0.60±0.02. 

PACS numbers: 64.70.Pf, 64.70.Dv, 71.15.Pd, 61.50.Ah, 61.66.f 



The use of computer models of liquids has played an 
important role for understanding glass forming-liquids. 
Computationally undemanding models are attractive, 
since long simulation times are important. Further more, 
the model does not have to be specific since the glass 
transition is universal. A simple model for molecular 
interactions is the the famous Lennard-Jones pair po- 
tential, Uij{rij) = 4.£{{rij/(JijY'^ - {^ij/^ijf)' It is not 
possible to investigate a single component Lennard-Jones 
liquid close to the glass transition since the structural re- 
laxation time is of the same order as the crystallization 
time. This can be avoided by using a binary Lennard- 
Jones liquid where the liquid consists of two kind of 
particles with differend radii cjij and binding energy £ij 
[U [2]. In 1991 Wahnstrom suggested parameters for a 
binary Lennard-Jones liquid in order to investigate the 
metastable undercooled liquid p. In the following years 
the model has been used by several others as a standard 

liquid lasiISllSlElElElIIo]. 

Increasing computer power expands the time scale that 
can be simulated. On a time scale of microseconds, the 
Wahnstrom liquid is no longer meta-stable, but crystal- 
lizes. In this paper we report observations of formations 
of crystallites and growth. It is possible to investigate 
a liquid, a highly viscous metastable liquid, a glass and 
crystallization within the same model. There is a sep- 
aration of phonons, structural relaxation, formation of 
crystallites, and crystal growth in time. 

The classical understanding of the crystallization pro- 
cess is as follows [11 . Crystallization consists of two 
events: formation of a critical crystallite and subsequent 
growth. The free energy of a crystallite (below the melt- 
ing temperature) consists of a positive surface part and 
a negative volume part. Crystallites below a critical size 
are unstable since the surface part will dominate. When 
a crystallite is above the critical size, the volume part 
will dominate, and the crystallite will grow. 

We have performed molecular dynamics simulations 
of the binary Lennard-Jones liquid suggested by Wahn- 



strom [T]. Particles interact via the Lennard-Jones pair 
potential. The system consists of two types of parti- 
cles labelled A and B. Energy is reported in units of 
£aa^ length in units of (Jaa^ mass in units of m^, tem- 
perature in units of T = Eaa/^Bi and time in units of 

= gaa \/ ^AA I ^AA ' Type B particles have a larger 
radius than type A particles, 0-33 = ^^aa and <jab = 
^{o'bb^(^aa) = l-lcTAA and twice the mass, tub = 2m^. 
All binding energies are equal, saa = £bb = ^ab- 

We define Argon units such that A particles may be 
considered to be Argon molecules; eaa = 1 kJ/mol, 
aAA 0.34 nm tua 40 g/mol, eAA/ks = 120.27 K 
and r = 2.15 ps. 

The simulations were performed using Gromacs soft- 
ware [12 . A switch function from 2.0(j aa to 3.5(Jaa was 
used to make potential go smoothly to zero at 3.5<jaa- 
Periodic boundaries were imposed. The Verlet velocity 
integrator with a time step of O.Olt* were used. The 
temperature was held constant using the Nose-Hoover 
thermostat [13 . The density were N/V = 0.75a^\ and 
Na — Nb = 512. This is the same density and Na/Nb 
ratio as originally used by Wahnstrom [1 , but twice the 
size. We also performed simulations with Na = '^Nb and 
a density of 0.82cr^^. The program VMD was used for 
visualization [M]. 

Fig.[l]shows the potential energy per particle as a func- 
tion of time for independent samples at T = 0.624 and 
T = 0.599. The composition was AB. At T = 0.624 the 
total simulation time covers 15 /is or in the order of IO^Tc^ 
structural relaxation times. One of the samples show a 
dramatic drop in energy due to crystallization. Crys- 
tallization is a rare event at this temperature but more 
common at the lower temperature. The samples that do 
not crystallize show some transient drops in energy. This 
can be explained as formation of small crystallites. This 
was confirmed by a common neighbor analysis (data not 
shown) [15j. 

The bottom panel on Fig. [l] show the quenched struc- 
ture of the crystallized sample. A crystal with a few 
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FIG. 1: The two top panels show the potential energy per 
particle as a function of time for independent runs at two 
temperatures with Na = Nb = 512. Several of the samples 
show a dramatic drop in energy due the growth of crystallites. 
At T — 0.624 the total simulation time of liquid is t = 13 fis 
(we made nine samples, not all data are shown), but only 
one sample crystallizes. The inset on the top panel shows a 
transient drop in energy which can be identified as temporary 
formation of a crystallite. At T = 0.599 crystallization hap- 
pens more often. The bottom panel the quenched structure 
of the crystallized sample at T = 0.624. Type A particles 
are green (dark gray) and type B particles are white (light 
gray). A crystal has formed with Na/Nb — 2, see Fig. [2j 
The excess B particles and some A particles form a liquid in 
the remained of the box. It was not possible to crystallize the 
whole sample within 4.5 /xs. 




FIG. 2: Part of crystal seen from different perspectives. The 
outline of a trigonal, hexagonal unit cell is sketched. Type A 
particles are colored in green (dark gray) and type B particles 
in white (light gray). The crystal structure is MgZn2 (C14) 
with P63/mmc symmetry. There are 6 type A particles on 
the h position, 2 on the a position and 4 type B particles on 
the / position (where /i, a and / are Wyckoff letters) [16 . 



defects has formed. It fills half the box, while the rest 
of the box is occupied by liquid. The crystal consists of 
more type A than type B particles so that Na/Nb — 2, 
whereas the liquid mainly consists of type B particles. 

Fig. [2] shows a part of the crystal seen from different 
perspectives. A unit cell is sketched. The crystal struc- 
ture is identified as the hexagonal MgZn2 Laves struc- 
ture. This is a close packing structure for the A2B com- 
position with (Jb/(^a = a/3/a/2 — 1.22 [T7 . This is in 
good agreement with the parameters used. 

Fig. [3] shows simulations of a perfect crystal in an or- 
thorhombic cell with 864 type A and 432 type B parti- 
cles. The density was p = 0.82. The crystal was heated 
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FIG. 3: The lower line shows the potential energy per particle 
during heating of a perfect crystal. The heating rates were 
2.2 X 10^ K//is and 2.2 x 10^ K//is in Argon units. At T = 1.3 
and T = 1.5 respectively the crystal melted. The systems 
were then cooled at the same rates. The liquids undergoes a 
glass transition at Tg ~ 0.6 (depending on cooling rate). The 
dashed line is a linear fit to the glass data with a slope of 
1.44. A partially melted crystal (approximately 90% crystal) 
was instantaneously cooled and kept at constant temperature 
T' for 30 ns. The x's mark the final energy at different T' 
temperatures. From this data, the melting temperature was 
estimated to be Tm — 0.93 =b 0.02. The dash-dot line have a 
slope of 3/2 corresponding to a classical harmonic crystal. 



at constant density from zero temperature with a rate 
of ^ = 4 x 10-4 and = 4 X IQ-^ At T = 1.3 
and T = 1.5 respectively the crystal melted. We expect 
the crystal to be over-heated at these temperatures. The 
melting temperature Tm was estimated by simulating a 
partially melted sample (approximately 90% of the vol- 
ume were crystal) at constant temperature T' for 30 ns. 
The sample melted when V > = 0.93 ± 0.02 and 
partially recrystallized when < T^. 

The liquid made from the melted crystal was cooled 
with a rate of ^ = -4 x 10"^ and ^ = -4 x lO"^. 
At Tm the sample did not crystallize, but became super- 
cooled. At Tg :^ 0.6 the sample passed though a glass 
transition (slightly dependent on cooling rate). 

In order to investigate the crystallization properties at 
the A2B composition further, we performed simulations 
of the super-cooled liquid at constant temperature with 
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FIG. 4: Potential energy per particle during simulations of 
metastable liquids at T = {0.625,0.583}, p = 0.82, Na = 
1024 and Nb = 512. Four independent samples are shown 
at each temperature. Six of the eight samples crystallizes. 
At T = 0.625 the limiting factor is the formation rate of 
crystallites, while growth rate is the limiting factor at T = 
0.583. 



N = 1536 and p = 0.82. Fig. |4] shows the potential 
energy per particle as a function of time at T = 0.624 
and T = 0.582. We used several independent samples. 
At the high temperature two out of four sample show a 
dramatic drop in energy. Here, the limiting factor for 
crystallization is the formation of a critical crystallite. 
When a critical crystallite have formed, the rest of the 
sample crystallizes fast. At the low temperature all four 
samples have formed a critical crystallite at t :^ but the 
growth rate is slow. From simulations at several temper- 
atures the maximum crystallization rate was estimated 
to T = 0.60 ± 0.02 (data not shown). 

In summary, we show for the first time that the Wahn- 
strom binary Lennard- Jones liquid undergoes crystalliza- 



tion when the metastable hquid is simulated for a time of 
lO^Ta close to the glass transition. The crystal structure 
is identified as MgZn2. This crystal is rather complex, 
but can form even when the composition is Na = Nb- 
The melting temperature of the crystal (at the A2B com- 
position and p = 0.82) is = 0.93 =b 0.02. The maxi- 
mum crystallization rate of the sample is T = 0.60 ±0.02. 
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